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1. Introduction 

The identification of the phases of QCD in the temperature - chemical potential plane and the 
exact localization and nature of the transitions among them are of central interest nowadays, due to 
the implications in cosmology, in astrophysics and in the phenomenology of heavy ion collisions. 
Most aspects of the QCD phase diagram are beyond the reach of perturbation theory and the lattice 
approach is a natural tool to face them. For non-zero chemical potential, however, the QCD fermion 
determinant becomes complex and the standard Monte Carlo importance sampling is unfeasible - 
the well-known "sign problem". 

Several strategies have been invented to circumvent this problem (for a review, see [|l|] and ^) 
and some useful information on the critical line separating the hadronic phase from the quark-gluon 
plasma phase in the region ii/T <\ has been already achieved. 

Here, we concentrate on one of these approaches, the method of analytical continuation, first 
used in Ref. ^ and in Ref. [Qj. The idea behind this method is very simple: numerical simula- 
tions are performed at imaginary chemical potential, Li = i/j.j, for which the fermion determinant 
is real, then Monte Carlo determinations are interpolated by a suitable function and finally this 
function is analytically continued to real values of Li. This method is rather powerful since j3 and 
li parameters can be varied independently and there is no limitation from increasing lattice size, 
as in methods based on reweighting. There is, however, an important drawback: the periodicity 
of the partition function and the presence of non-analyticities arising for imaginary values of the 
chemical potential [§] make that the region useful for numerical determinations is inside the strip 
< jJ-i/T < n/3. This implies that the accuracy in the interpolation of the results at imaginary 
chemical potential has a strong impact on the extension of the domain of real /I values reachable 
after analytical continuation. 

So far, the method of analytical continuation has been applied in SU(3) with rif = 2 [§J, rif = 
3 [Q] and rif = 4 [||]. Moreover, it has been tested in several theories which do not suffer the sign 
problem, by direct comparison of the analytical continuation with Monte Carlo results obtained at 
real ii 

In all these applications a truncated Taylor series (or, more simply, a polynomial) has been 
used as interpolating function. In this work we show that the method can be considerably improved 
if ratios of polynomials are used instead. As a test-field for our proposal we adopt the same theory 
used in Ref. [Q], namely SU(2) (or 2-color QCD), with rif = 8 degenerate staggered fermions. 

2. Theoretical background 

Long ago Roberge and Weiss have shown [j|] that the partition function of any SU(A0 gauge 
theory with non-zero temperature and imaginary chemical potential, /I = iili, is periodic in 6 = 
Lii/T with period 2n/N and that the free energy F is a regular function of 6 for T < Te, while it is 
discontinuous at 6 = 2%{k + 1 /2)/N, k = 0, 1,2, . . ., for T > Te, where Te is a characteristic tem- 
perature, depending on the theory. The resulting phase diagram in the (T, 0)-plane is given in Fig. |T] 
(left), where the vertical lines represent first order transition lines. This structure is compatible with 
the li — ► —Li symmetry, related with CP invariance, and with the Roberge-Weiss periodicity. The 
/^/-dependence of any observable is completely determined if this observable is known in the strip 
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Figure 1: (Left) Phase diagram in the (T. 9) plane according to Ref. [Q]. (Right) Tentative phase diagram in 
the (T, 9) plane after the inclusion of the chiral critical lines. 

< 6 < 71 /N. These predictions have been confirmed numerically in several cases, studying the 
behaviour of quantities like the Polyakov loop and the chiral condensate |E| §, §]. 

A phase diagram like that in Fig. [I] (left) would imply the absence of any transition along the 
T axis in the physical regime of zero chemical potential for any value of N, of rif and of the quark 
masses, which cannot be true. Therefore, it is necessary to admit that the phase diagram in the 
(r, 0)-plane is more complicated than in Fig. [T] (left). The simplest possibility is given in Fig. |T] 
(right), where the added lines generally represent transitions which can be first order, second order 
or crossover. The temperature T c is the critical or pseudo-critical one for the transition at zero 
chemical potential. It is convenient to redraw the phase diagram of Fig. [j] (right) in the (j8,/t/)- 
plane (Fig. |J), where j8 = 2N/g 2 , /t/ = a\ii is the imaginary chemical potential in lattice units and 
it has been used the fact that T = 1 /(aN T ), with N T the temporal extension of the lattice. 

3. Numerical results 

We performed numerical simulations on a 16 3 x 4 lattice of the SU(2) gauge theory with 
rif = 8 degenerate staggered fermions having mass am = 0.07. For this theory the tentative phase 
diagram looks like in Fig. §, with f} E — 1.55 |§] and j8 c ~ 1.41 [01- We used the hybrid Monte 
Carlo algorithm, with dt = 0.01. The observables we determined are the Polyakov loop, the chiral 
condensate and the fermionic number density. We have chosen j8 = 1.90 > (5e and have taken 
1000-5000 measurements for several values of the imaginary chemical potential in the interval 
from zero to the value corresponding to the first RW transition line, fii = 7l/8, and for several 
values of the real chemical potential ji. R in the interval [0, 2] . Simulations have been performed on 
the APEmille crate in Bari and on the recently installed computer facilities at the INFN apeNEXT 
Computing Center in Rome. 

We have used the data at imaginary chemical potential jij to determine the parameters of the 
interpolating function, then we have analytically continued this function to real values of the chem- 
ical potential and compared there with direct Monte Carlo determinations. In order to fulfill CP 
invariance, the interpolating function must be a function of ji 2 for observables, such as the Polyakov 
loop and the chiral condensate, which do not depend explicitly on /x. The fermion number density, 
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Figure 2: Phase diagram in the (/3 , /i/)-plane; N is the number of colors, N r the extension of the lattice in the 
temporal direction. The numerical values for Pe and j3 6 are valid for SU(2) in presence of ttf — 8 degenerate 
staggered fermions with mass am — 0.07. 



being the logarithmic derivative of the partition function with respect to the chemical potential, is 
instead an odd function of LL. For the Polyakov loop and the chiral condensate we have considered 
a second order polynomial in LL 2 , 

A + BLL 2 + Cp 4 , (3.1) 
according to the standard approach, and the ratio of two first oder polynomials in ll 2 , 

A + Bu 2 

TW ' (3 " 2) 

according to our new proposal. Similarly, for the fermionic number density we have used a poly- 
nomial of the form 

Ah+Bli 3 +Ch 5 , (3.3) 

and the ratio 

Our findings are summarized in Figs. ^ ^ and 0. In Fig. || we put on the same plot the imag- 
inary part of the fermionic number density as a function of LLj and the real part of the fermionic 
number density as a function of the real ll. The two data sets match smoothly in LL = 0, which is 
a necessary condition for the applicability of the method of analytical continuation. The fermionic 
density approaches 2 for large values of the real chemical potential. This saturation effect is arti- 
ficial and is due to the fact that no more than two fermions per site can be accommodated in the 
lattice ("Pauli blocking"). The solid lines represent the two kinds of interpolating functions, whose 
parameters are determined by a fit on the data at imaginary chemical potential. Here both interpo- 
lations, polynomial and rational function, nicely reproduce the data at real ll over a large interval. 
Deviations start at values of \i for which the saturation effects are certainly important. 

In Fig. ^| we show the chiral condensate as a function of ll 2 . Again data at imaginary LL, i.e. 
LL 2 < 0, and data at real }x, i.e. LL 2 > 0, nicely match at /i = 0. This time the different behavior of the 
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Figure 3: Negative side of the horizontal axis: imaginary part of the fermionic number density vs. the 
imaginary chemical potential. Positive side of the horizontal axis: real part of the fermionic number density 
vs. the real chemical potential. The green (blue) solid lines represent the polynomial (ratio of polynomials) 
interpolating function; the dashed lines give the corresponding uncertainty, coming from the errors in the 
parameters of the fit. 



two kinds of interpolation clearly emerges. The ratio of first order polynomials in \i 2 reproduces 
the data at real /I on a much larger interval than the second order polynomial in \i 2 . Deviations 
arise for values of real /I for which saturation effects are probably already important. The same 
conclusions can be drawn from Fig. |5] which shows data and interpolations for the Polyakov loop. 

The above conclusions do not change if larger order terms are included in the polynomial 
interpolation (3J_). In fact, larger order polynomials fail to reproduce the data at real \i even earlier 
in \i than second order polynomials. This is due to the fact that the higher order terms of the 
polynomial are the less accurately determined in the fit to data at imaginary \i. On the other side, if 
in the ratio of polynomials the order of the polynomials at the numerator and/or at the denominator 
is increased, no improvement is observed. 

An approach alternative to that based on ratio of polynomials is to use large order polynomials 
and then to build out of them rational functions by means of Pade approximants 1 . This alternative 



1 The use of Pade approximants for the analytical continuation of the critical line from imaginary to real chemical 
potential has been suggested in Ref. [112]. 
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Figure 4: Chiral condensate vs. \xr. The solid and dashed lines have the same meaning as in Fig. 3. 



approach behaves in a way comparable to ratio of polynomials. 

4. Conclusions and outlook 

We have verified by comparison with direct Monte Carlo determinations at real chemical po- 
tential in 2-color QCD that the method of analytical continuation considerably improves if ratio of 
polynomials is used as interpolating function instead of truncated Taylor series. 

In the case of the Polyakov loop and of the chiral condensate an interpolation of numerical 
data at imaginary chemical potential over the window permitted by Roberge-Weiss singularities 
allows an extrapolation to real values of the chemical potential over a much larger region. 

Deviations at very large values of the chemical potential could be due to unphysical saturation 
of the fermionic density ("Pauli blocking"). 
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